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ABSTRACT 

We report the detection of radio emission from PSR J1311— 3430, the first milhsecond pulsar dis- 
covered in a bhnd search of Fermi Large Area Telescope (LAT) gamma-ray data. We detected radio 
pulsations at 2 GHz, visible for <10% of ^^4.5 hours of observations using the Green Bank Telescope 
(GBT). Observations at 5 GHz with the GBT and at several lower frequencies with Parkes, Nangay, 
and the Giant Metrewave Radio Telescope resulted in non-detections. We also report the faint de- 
tection of a steep spectrum continuum radio source (0.1 mJy at 5 GHz) in interferomctric imaging 
observations with the Jansky Very Large Array. These detections demonstrate that PSR J1311— 3430 
is not radio quiet, and provides additional evidence that radio-quiet millisecond pulsars are rare. The 
radio dispersion measure of 37.8 pc cm~'^ provides a distance estimate of 1.4 kpc for the system, 
yielding a gamma-ray efficiency of 30%, typical of LAT-detected MSPs. We see apparent excess delay 
in the radio pulses as the pulsar appears from eclipse and we speculate on possible mechanisms for 
the non-detections of the pulse at other orbital phases and observing frequencies. 
Subject headings: pulsars: individual (PSR J1311— 3430) 



1. INTRODUCTION 

The Large Area Telescope (LAT) on the Fermi 
Gamma-ray Space Telescope has been surveying the 
gamma-ray sky in the 20 MeV to 300 GeV band since 
August 2008. The LAT survey has vastly increased the 
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population of known gamma-ray emitting sources includ- 
ing pulsars, blazars, supernova remnants, pulsar wind 
nebulae, and more. Howev er, in the second LAT source 
catalog (INolan et al.l[2nTl hereafter 2FGL) 575 (31%) 
of the 1873 sources were not associated with any known 
counterpart. The properties of the unassociated source 
population have been studied (jAckermann et ahl l2012f) 
and they roughly fall into two categories of 'pulsar-like' 
(non- variable with spectra that exhibit significant curva- 
ture) and 'blazar-like' (highly variable, spectra without 
strong cutoffs). 

Beyond studying their gamma-ray emission properties, 
the most promising way to make progress on understand- 
ing and identifying the unassociated sources is through 
multi-wavelength observations. As a striking example, 
deep searches for radio pulsars in the pulsar-like popu- 
lation have led to the discovery of 43 new millisecond 
pulsars (MSPs) and 4 young and middle-aged pulsars 
(|Rav et all 120121) . Of course, since all of these pulsars 
were discovered in radio searches, these discoveries do not 
probe the possible radio-quiet MSP population, which 
must be discovered in other wavelengths. 

The first strong candidate radio-quiet MSP was found 
through X-ray and optical studies of the bright unassoci- 
ated source 2FGL J2339. 6-0532 l|Romam fc Shawll2011l : 
Kong ct al. 2012). The orbital modulation of the optical 
emission sugge sts the system comprises a "black widow" - 
type MSP (King et al.l l2005[ ) in an 4.6-hr orbit around 
a low-mass companion. Very re cently, radio puls ations 
from this source were reported (|Rav et al.l I2012D . con- 
firming that it is indeed an MSP whose wind is evapo- 
rating its compa nion. 

More recently, iRomanil (|2012[ ) identified another such 
system associa ted with 2FGL J1311 .7— 3429 (hereafter 
J1311; see also iKataoka et all (|2012l )). In this case, the 
optical modulation revealed an orbital period of 1.56 
hours, the shortest of any known MSP, with a compan- 
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ion that is very hydrogen poor. Using the optical source 
position and information abou t the orbital parani eters 
from the optical measurements. I Pletsch et alj (|2012[ ) dis- 
covered the first gamma-ray MSP found in a blind search, 
PSR J1311-3430, with period P = 2.56 ms. This dis- 
covery confirmed the black widow pulsar nature of the 
system and yielded important parameters such as the 
spin-down luminosity {E = 4.9 x 10"^* erg s~^) of the 
pulsar and the precise orbital parameters, including the 
projected semi-major axis of the orbit (10.6 It-ms). As 
the pulsar was discovered in the gamma-ray band, it was 
a candidate for being the first radio-quiet MSP. Analy- 
sis of archival 820 MHz and 350 MHz radio observations 
with the Robert C. Byrd Green Bank Telesc ope (GET) 
resulted in no det ections (jPletsch et al.lf20T^ . However, 
as Romani! ()2012t ) pointed out, the low frequency radio 
emission could be scattered or absorbed and thus be un- 
detectable from Earth, even if the pulsar radio beam was 
directed at Earth. Because scattering and absorption of 
the radio emission are strong functions of frequency, this 
motivates further radio observations at higher frequen- 
cies. 

In this Letter, we report the detection of radio emis- 
sion from PSR J1311— 3430, both in continuum imaging 
observations with the Karl G. Jansky Very Large Array 
(VLA) and in a radio pulsation search with the GBT at 
2 GHz. This confirms that PSR J1311-3430 is indeed a 
radio MSP with the beam sweeping across Earth, com- 
pleting the chain from unassociated gamma-ray source 
to optical identification to gamma-ray pulsation detec- 
tion to radio detection. 

2. OBSERVATIONS AND RESULTS 

2.1. VLA Imaging Observations 

We observe d the J1311 fie ld on 2011 January 14 with 
the VLA (Perl ev et all 1201 If ) while in the C -array (pro- 
gram S2270). Observations were made at 1.4 and 5 GHz 
using two 128 MHz wide intermediate frequencies cen- 
tered at 1.327 + 1.455 GHz and 4.895 -t- 5.023 GHz, 
respectively. In the 1-hr observing run, we obtained two 
9-min target scans per frequency bracketed with scans of 
a phase calibrator (J1316— 3338), while the primary flux 
calibrator was 3C 286. 

An analysis of the 5 GHz data revealed an unre- 
sol yed (rat i o of t he peak to integrated flux density ~ 
1; iCondonI (|1997[ )) 0.10 ± 0.04 mJy point source at a 
position (J2000.0), R.A. = 13'>ll'°45^78 (±1.0") and 
Deck = -34°30'31.2" (±1.7"), consistent with the pul- 
sar (i.e. the pulsar is within the VLA synthesized beam). 
In the 1.4 GHz image, there is a 0.99 ± 0.44 mJy source 
(peak = 0.33 ± 0.11 mJy/beam) that shows evidence of 
extension, with a centroid that is slightly offset (~ 13") 
from the 5 GHz position (Figure [ij . Imaging the two 
scans of J1311 at each frequency separately (see Figure[2] 
for the corresponding orbital phases) did not reveal any 
statistically significant fiux changes. 

2.2. Radio Pulsar Observations 

After the discove ry of the gamma-ray pulsations by 
IPletsch et al.l ()2012l ). we performed a series of radio ob- 
servations to search for radio pulsations from the pul- 
sar. Given the non -detections at 820 and 350 MHz 
(jRansom et al.l 120111: [Pletsch et al..,2012.) . and the 1.4 
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Fig. 1. — VLA radio contours at 5 GHz (white; 80 and 113 
/ijy/beam) and 1.4 GHz (green; 250 and 354 /ijy/beam) over- 
laid onto a Chandra 0.5 — 7 keV image of the region showing the 
X-ray point source, CXOU J131145. 7-343030 from Cheung et all 
II2012I ). found coincident with the pulsar. The synthesized beams 
are 16.2" x3. 5" (position angle, PA = -26.1 deg) and 49.1" X 14.7" 
{PA = —22.9 deg), respectively. 
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Fig. 2. — Figure showing the orbit of the pulsar with the phase 
range of the 2 GHz radio pulsation detection marked in gray. The 
VLA observation intervals centered at phases 0.51 and 0.87 (5 GHz; 
cyan) and 0.64 and 1.00 (1.4 GHz; red) are marked. Orbital phase 
is measured as a fraction of the orbit from the ascending node 
(when the pulsar crosses the plane of the sky heading away) ; orbital 
phases of 0.0, 0.25, 0.5, and 0.75 are noted. Superior conjunction 
(e.g. mid-eclipse time) is also shown. 



and 5 GHz VLA detections (t j2.ip . the priority was on 
higher frequency observations to minimize the effects of 
scattering, dispersion, and absorption. For completeness, 
we include the earlier GBT observations in our analysis. 
The observations performed are listed in Table [T] All ob- 
servations are made with 2 summed polarizations. In all 
cases, we compute the pulsed flux density limit s using the 
modifi ed radiometer equation as described bv lRav et al.l 
(|2011f ). We assume a pulse duty cycle of 0.1 and a signal 
to noise threshold for detection of 5a. 
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2.2.1. Green Bank Telescope 

We observed the J1311 field with the GBT (project 
GBT/12A-487) with the GUPPf backend (DuPlaiireril] 
[2008. ) on 2012 July 28 for 1.4 hours at 2 GHz and 1.3 
hours at 5 GHz with a bandwidth of 800 MHz and 40.96 
/iS sampling. At 2 GHz we used 2048 frequency channels, 
while at 5 GHz we used 1024 channels. We obtained a 
second GBT observation on 2012 August 19 spanning 
3 hours at 2 GHz, with the same observing setup. At 
the start of each observation, we observed the pulsar 
B1257+12, which was detected each time^ 

We analyzed the data using PRESTO (jRanso ^ [200lfl 
to first excise strong radio frequency interference (RFI) 
signals and then fold the data using an ephemeris based 
on the parameters of the gamma-ray pulsation discov- 
ery. The RFI environment during these observations was 
rather benign, requiring blanking of some narrow band 
signals amounting to < 7% of the data. The only free 
parameter in the search was dispersion measure (DM). 
In our first 2 GHz observation, we detected radio pul- 
sations during an ^ 1100-s interval beginning at MJD 
(UTC) 56136.908 with a DM of 37.84 ± 0.26 pc cm'^ 
(see Figure [3]). The significance of the detection over the 
1100 seconds where the pulse is visible, after correcting 
for the 6400 trials in the DM search, is 10.3 cr. This signif- 
icance is computed from the profile using the test for 
excess variance. It differs from the value shown in Figure 
[3] because it is based on a subset of the data that exclude 
the first 600 seconds where no signal is apparent. Based 
on the detection significance and the sensitivity of our 
observations, we estimate a mean flux density of 60 ± 30 
/ijy. No pulsations were seen in the 5 GHz observation 
or the second 2 GHz observation. 

As seen in Figure [3l the signal significance peaks at 
a frequency slightly offset {Av ~ 5 x 10"'^ Hz) from 
that predicted by the gamma-ray ephemeris. The orbital 
phase at the start of the interval in Figure[3]is 0.24, when 
the pulsar is expected to be at mid-eclipse, so this is ap- 
parently caused by excess delay in the radio pulse as the 
pulsar is emerging from eclipse, as is seen in many eclips- 
ing MSPs. To test this, we folded only the data starting 
after 600 seconds and the detection significance peaks at 
the predicted period and period derivative. In addition, 
we generated 15 times of arrival (TOAs) across the ob- 
servation, with 114 seconds of integration per TOA. We 
discarded the first three with no apparent pulsed detec- 
tion. Looking at only the last 10 TOAs, the measured 
period is consistent with that predicted from the gamma- 
ray ephemeris. The two earlier TOAs appear to show a 
changing delay that is about 0.6 ms (1/4 of a pulse pe- 
riod) just as the pulse becomes visible. 

2.2.2. Parkes 

We made a total of 7 observations at 1.4 GHz us- 
ing the Parkes Telescope with the Analog Filter Bank 
(AFB) backend recording 256 MHz of bandwidth (0.5 
MHz channels sampled every 125 /is) and analyzed them 
in the same manner as described above. No pulsations 
were detected in any of the observations. 

2.2.3. Nangay 

We performed 5 observations with the Nangay radio 
telescope (NRT) around 1.5 GHz ranging from 700 to 



TABLE 1 

Radio Pulsar Observations 



Obs Code 


Date 


^obs 


Orb Phase 




c b 
'-'mill 








(s) 






(K) 


(^Jy) 


nRT-S2n 


2011 


Jan 15 


2577 


0.092 - 


- 0.550 


38 


53 (22) 


nRT-S2n 


2011 


Dec 10 


1723 


0.346 - 


- 0.652 


38 


53 (22) 


nRT ^^90 

V_JJ_> ± C>^\J 


2012 


Aug 24 


900 


0.295 - 


- 0.455 


38 






2012 


Aug 25 


900 


0.977 - 


- 0.137 


38 






2012 


Feb 24 


1287 


0.339 - 


- 0.568 


104 


203 (22^ 


GBT-S 


2012 


Jul 28 


4939 


0.859 - 


- 0.737 


26 


114 it 60 


GBT-S 


2012 


Aug 19 


11697 


2 orbits 


26 


20 ('36~1 


GBT-S 


2012 


Aug 25 


900 


0.346 - 


- 0.506 


26 


20 (36^ 


GBT-C 


2012 


Jul 28 


4552 


0.931 - 


- 0.740 


22 


16 (115) 


GMRT-fifiy 


2011 


Feb 15 


3386 


0.108 - 


- 0.710 


108 


445 (117) 


GMRT-BO? 


2011 


Jun 23 


7246 


0.349 - 


- 0.637 


108 


445 (117) 


GMRT-322 


2012 


Aug 15 


3607 


0.844 - 


- 0.485 


137 


184 (18) 


Parkes- AFB 


2012 


Mar 22 


3600 


0.606 - 


- 0.246 


29 


115 (114) 


Parkes- AFB 


2012 


Jul 11 


3600 


0.337 - 


- 0.977 


29 


115 (114) 


Parkes- AFB 


2012 


Jul 12 


3600 


0.025 - 


- 0.665 


29 


115 (114) 


Parkes- AFB 


2012 


Jul 16 


3600 


0.092 - 


- 0.732 


29 


115 (114) 


Parkes- AFB 


2012 


Aug 17 


1635 


0.403 - 


- 0.694 


29 


115 (114) 


Parkes- AFB 


2012 


Aug 17 


5357 


0.766 - 


- 0.718 


29 


115 (114) 


Parkes- AFB 


2012 


Aug 24 


4500 


0.130 - 


- 0.930 


29 


115 (114) 


NRT-L 


2012 


Jul 16 


2512 


0.182- 


-0.628 


39 


48 (54) 


NRT-L 


2012 


Jul 21 


2997 


0.563- 


-0.096 


39 


48 (54) 


NRT-L 


2012 


Aug 23 


736 


0.368- 


-0.499 


39 


48 (54) 


NRT-L 


2012 


Sep 1 


2996 


0.788- 


-0.320 


39 


48 (54) 


NRT-L 


2012 


Sep 24 


1257 


0.094- 


-0.317 


39 


48 (54) 



^ System temperature including receiver temperature (Tree) and 
sky temperature (T^j^y) from scaling the 408 MHz Haslam map to 
the observing frequency and adding 2.7 K for the cold sky back- 
ground. 

^ Flux density limits are computed for a standard integration time 
of 1000 s at the observing frequency (corresponding to the duration 
where pulses were seen in the one detection), with the equivalent 
limit at 1400 MHz in parentheses, assuming a spectral index of 
— 1.8. The bold value is the detected flux density in this observa- 
tion, scaled to 1.4 GHz. 

3000 s. We used a 5 12- MHz bandwidth divided into 
1024 channels with 64-/is sampling. Dedispersion at a 
DM of 37.84 pc cm"'^ and folding using the known pa- 
rameters did not produce any detectable pulsations. The 
offsets from the gamma-ray pulsation position are negli- 
gible when compared to the beam size of the telescope. 

2.2.4. Giant Metrewave Radio Telescope 

We analyzed two archival 607-MHz (32 MHz band- 
width) observations made using the Giant Metrewave 
Radio Telescope (GMRT) incoherent array in early 2011 
as part of a pulsar search of LAT unassociated sources. 
Folding these data with the timing model from the LAT 
discovery resulted in no detections. Later with the pre- 
cise pulsar position, we made sensitive coherent array 
observations at 3 22 MHz, also with 32 MHz bandwidth 
(|R,ov et al.ll20T(il) . A gain, no pulsations were detected. 

2.3. Radio /Gamma- Ray Profiles 

We computed phase-aligned radio and gamma-ray pro- 
files, shown in Figure |4l The gamma-ray profiles are 
comp uted from the same weighted LAT events as in 
iPletsch et al.l ()2012[ ) spanning 2008 August 4 to 2012 
July 10. The ph ase was compu t ed ac cording to the pul- 
sar ephemeris of IPletsch et al.l (|2012[ ). with the fiducial 
phase shifted to correctly align with the radio profile, 
using a dispersion measure of 37.84 pc cm~^. The radio 
profile at 2 GHz was constructed from the GBT observa- 
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2 Pulses of Best Profile 



^J2000 ■ 



Search Information 
RAj2(j(,o = 13:11:45.7440 DEC 

Folding Parameters 
Reduced = 4.696 P(NoIse) < 9.03e-17 (»,8.2ct) 
Dispersion Measure (DM; pc/cm"^) = 37.84 



34:30:30.2400 



2.56059866(20) 
0.0(9. 2)x10"'- 
= 0.0(3.5)x10" 
Binary Parameters 
Porb = N/A 
a^sln(l)/c (s) = N/A 



bary 

P'V... (s/s^ 



N/A 
= N/A 
= N/A 




Phase 
5136_131 1- 



Reduced x' 
■3430_0002_0011.flti 



(pc/cm ) 



2.56059866 (ms) 

pauir 19-0ct-2012 



Fi g. 3. — Radio detecti on of PSR J1311— 3430 using the GBT at 2 GHz. This shows a span of 1711 seconds folded using the ephemeris 
from iPTetsch et al.l 1)20121 ). The pulsar is very weak during the first 600 seconds of this interval and appears to show some excess delay as 
it comes out of eclipse. When the first 600 seconds are excluded from the analysis, the detection significance is above 10 tr. The fact that 
the detection is strongest in a frequency band of ~ 150 MHz and weak in the rest of the band is consistent with interstellar scintillation. 



tion at 2 GHz using the 1100 seconds of data where the 
pulsar signal is apparent. 

Both the radio and gamma-ray profiles show a double 
peak structure and are roughly aligned. In comparing the 
radio to gamma-ray profile alignment, it is important to 
consider the uncertainty in absolute phasing due to the 
uncertainty in the DM measurement. The uncertainty 
in the DM is 0.26 pc cm"'^. A change in the DM of this 
magnitude corresponds to a shift in the relative phasing 
of the radio and gamma-ray profiles of 0.1 in phase. So, 
within the uncertainties, either the first radio peak could 
be aligned with the first gamma-ray peak, or the second 
peaks could be aligned. However, the peak separation in 
the radio profile is 0.28 ± 0.02 and the gamma-ray peaks 
are separated by 0.35 ± 0.01, so it not possible for both 
peaks to be precisely aligned. Finally, the radio profile 
also may be broadened by scattering but without a higher 
quality profile, or detections at other frequencies, this is 
hard to quantify. 

3. DISCUSSION 

This detection demonstrates that although this pul- 
sar would have been extraordinarily difficult to detect 
in a radio search, it is indeed visible from Earth in the 
radio band, at least when the conditions in the system 
permit. The measurement of the pulsar dispersion mea- 
sure yields a distance estimate fr om the NE2001 galac- 
tic free electron density model (jCordes &: Laziol |2002|) 



of 1.4 kpc, typical of the LAT-detected MSP sample. 
Using this distance and the measured proper motion 
of 8 mas yr~^ we compute that the contribution to 
the observed per iod derivative from the Shklovskii effect 
(jShklovskiil 119701 ) is only 2.7%. Using the DM distance 
to compute the gamma-ray luminosity gives an efficiency 
rj = Lry / E = 30%, assuming that the beaming correction 
fn is 1. 

To gauge PSR J1311— 3430's signi ficance f o r the 
gamma-ray pulsar population, we follow iRomanI (poTl ) 
and restrict our attention to the brightest ~ 250 Fermi 
sources, for which extensive counterpart studies (in- 
cluding repeated radio pulse searches) have been made. 
These include 14 gamma-ray MSPs, all of which have 
now been detected in radio plus the newly-discovered 
PSR J2339— 0533. There are four remaining unassoci- 
ated objects in this sample - their gamma-ray proper- 
ties (a significant high-energy spectral cut-off and low 
spectral variability) mark them as likely pulsars. How- 
ever some are at lower Galactic latitudes and so, unlike 
J1311— 3430 and J2339— 0533, may be powered by young 
pulsars. Even if we count all these sources as MSPs un- 
detectable in the radio, the fraction is only 4/18 = 0.22. 
This is dramatically lower than the 27/42 = 0.64 frac- 
tion of radio-quiet young pulsars in this sample. We con- 
clude that, for MSPs, radio beams cover nearly as much 
sky as the gamma-ray beams, and largely overlap. This 
agrees with previous results showing that the beaming 
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Fig. 4. — Phase-aligned radio and gamma-ray profiles. The 
gamma-ray profiles are weighted counts with energies > 100 MeV 
and the radio profile is from the GBT detection of pulsations at 2 
GHz. The dotted lines are the same profile phase shifted forward 
and back by 0.1 in phase, comparable to the DM uncertainty. 

fraction of MSP s is near unity and likely from a fan beam 
(|Lorimed I2008L and references therein). Conversely the 
lack of a large number of black- widow-like LAT sources 
remaining to be detected suggests that there is not a large 
phase space for binary gamma-ray pulsars that cannot 
(at least occasionally) be detected in the radio; by the 
time the wind from the companion is strong enough to 
completely bury the radio signal, the pulsar accelerator 
itself may be quenched. 

The small fraction of time that the radio pulsations are 
detectable is notable. There are several possible explana- 
tions that are difficult to differentiate between with only 
one relatively low signal-to-noise detection. One reason 
could be simply interstellar scintillation. High frequency 
observations of nearby MSPs are particularly susceptible 
to scintillation. For some of the MSPs found in LAT- 
directed surveys at Parkes, the observed flux density 
varies by a factor of 20 from observation to observation 
and some of the pulsars are detectable in only about one 
third of the observations (Camilo et al. 2012, in prepara- 
tion). Alternatively, the pulsations could be eclipsed by 
material local to the system itself. Ther e are several pos- 
sible mechanisms for such eclipses (see [Thompson et al.l 
11994, for a review) and which mechanism is operative 
may be different for different systems. In this case, the 
continuum detections of a radio point source suggest that 
the pulsations may be scattered out of existence (e.g. 
scattering by plasma turbulence) rather than the radio 
waves being actually absorbed. Variations in the scat- 
tering medium may be responsible for the non-detections 



at other times and frequencies. The presence of a red, 
flaring optical component, likely fro m the pulsar flux re- 
processed in the companion wind (i Romanf "eFalllMl) 
gives additional evidence for a wind of variable density 
and covering fraction. The timescale and magnitude of 
dispersion measure variations can be a useful diagnostic 
of these processes. We see some evidence for DM grow- 
ing to just over 38.0 pc cm~'^ in the later part of our 
observation but the significance is too limited to draw 
any firm conclusions. Also, it is difficult to constrain the 
presence of a scattering tail in the pulse profile for the 
same reasons. 

Considering the VLA imaging observations, if a sig- 
nificant portion of the flux at 1.4 GHz is from the point 
source, it would imply a spectral index of —1.8, typical of 
known pulsars (Maron et al. 2000) , although there is an 
uncertain contribution from, e.g., a bow shock or pulsar 
wind nebula. Deeper imaging in a higher resolution array 
configuration is required to characterize the continuum 
source and its variability properties. 

Within the uncertainties in the absolute phasing, the 
radio and gamma-ray pulse profiles are approximately 
aligned, although their peak separations are slightly dif- 
ferent, with the radio peaks being closer together than 
the gamma-ray ones. This could indicate that both are 
formed in the same geometric region, but at different 
altitudes. 

This detection illustrates that additional observations, 
particularly at high frequency, are needed to search for 
radio pulsations from the other strong radio-quiet MSP 
candidates. Of course, these searches would be much 
more sensitive with a known pulsation from a gamma- 
ray detection, as was the case with J1311. Exhaustive 
searches of the LAT data as well as repeated radio ob- 
servations are clearly warranted in these cases. 
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